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Abstract —In this paper, a new cooperation structure for 
spectrum sensing in cognitive radio networks is proposed which 
outperforms the existing commonly-used ones in terms of energy 
efficiency. The efficiency is achieved in the proposed design by 
introducing random interruptions in the cooperation process 
between the sensing nodes and the fusion center, along with 
a compensation process at the fusion center. Regarding the 
hypothesis testing problem concerned, first, the proposed system 
behavior is thoroughly analyzed and its associated likelihood-ratio 
test (LRT) is provided. Next, based on a general linear fusion rule, 
statistics of the global test summary are derived and the sensing 
quality is characterized in terms of the probability of false alarm 
and the probability of detection. Then, optimization of the overall 
detection performance is formulated according to the Neyman- 
Pearson criterion (NPC) and it is discussed that the optimization 
required is indeed a decision-making process with uncertainty 
which incurs prohibitive computational complexity. The NPC is 
then modified to achieve a good affordable solution by using 
semidefinite programming (SDP) techniques and it is shown that 
this new solution is nearly optimal according to the deflection 
criterion. Finally, effectiveness of the proposed architecture and 
its associated SDP are demonstrated by simulation results. 

Index Terms —Cognitive radio (CR), cooperative spectrum 
sensing, efficiency, decision/data fusion, correlation, non-ideal 
reporting channels. 


I. Introduction 

E FFICIENT utilization of available resources is a vital 
requirement in designing modern communication sys¬ 
tems. In order to increase the efficiency, significant levels of 
flexibility and adaptability are realized in wireless networks 
by employing the cognitive radio (CR) technology which is 
developed through adding certain artificial-intelligence-based 
capabilities, such as self-awareness, context-awareness, and 
machine learning to software-defined radios |1J. Due to these 
prominent capabilities, CR has been of great interest for de¬ 
veloping green and spectrum-efficient wireless communication 
systems, see e.g., |2| and references therein. 

Reliable spectrum sensing is of major importance in de¬ 
signing CR networks (CRNs). The average throughput of the 
secondary users (SUs), their consumed energy, and the amount 
of interference experienced by the primary users (PUs) are 
directly related to the effectiveness of the sensing methods 
incorporated in CRs. In fact, a reliable sensing algorithm 
is necessary to obtain suitable transmission opportunities 
without compromising the integrity of the primary network, 
particularly in low signal-to-noise ratio (SNR) regimes |3j. 
Nevertheless, impairments such as shadowing and multipath 
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fading associated with typical wireless environments often 
degrade the spectrum sensing performance. In order to mitigate 
this issue, the reliability of the sensing process is significantly 
enhanced through establishing cooperation among spatially- 
divers sensing nodes. This cooperation, referred to as coop¬ 
erative spectrum sensing, can be realized in different config¬ 
urations among which the so-called centralized cooperative 
sensing is a commonly-used one J2J. In this structure, the 
cooperation is coordinated by and the overall sensing outcome 
is generated in a special node called the fusion center (FC), see 
e.g.. This cooperation is in fact performed as a three- 

phase process, i.e., each node first performs spectrum sensing 
individually by using its own built-in sensing scheme. Then, 
the sensing nodes send their local sensing outcomes to the FC 
through the so-called reporting channels and finally, the FC 
combines the received local sensing outcomes to decide the 
presence or absence of the PU. 

Censoring m-m is a common design strategy for reduc¬ 
ing the energy consumption of cooperative spectrum sens¬ 
ing schemes. In this method sensors are assumed to censor 
their observations so that each sensor sends to the FC only 
informative observations, and leaves those deemed uninfor¬ 
mative untransmitted. In hard-decision-based cooperation with 
censoring 1131—1161, instead of comparing the local sensing 
outcome with a single threshold to make the binary decisions, 
it is compared against two thresholds, i.e., an upper threshold 
and a lower one. If the sensing outcome is below the lower 
threshold, the sensing node reports the null hypothesis by 
sending, let’s say, 0 to the FC and if the local sensing outcome 
is above the upper threshold, then the sensing node reports 
the alternative hypothesis by sending 1 to the FC. However, if 
the sensing outcome resides between the two thresholds, it is 
not considered reliable enough to make a local decision and 
the CR node does not report anything to the FC, saving the 
energy and communication resources. In soft decision based 
cooperation 0-GiJ each sensing CR node compares the 
sensing outcome with a predefined threshold and if it is below 
the threshold, then the sensor avoids reporting to the FC. 

Sensor selection and sensor scheduling | 211 © 

algorithms provide a different set of methods for improving 
the energy efficiency in distributed detection schemes. Sensor 
selection algorithms mainly aim at selecting a set of sensors 
whose collaboration contributes significantly to the overall 
detection process. The goal of sensor scheduling algorithms 
in general is to increase the energy efficiency of cooperative 
sensing by dividing the sensors into non-disjoint groups which 
activate successively. In other words, only the sensors from the 
current active set are responsible for sensing and for reporting 
the sensing outcome while all other nodes are in a low-energy 
sleep mode. 
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A. Related Work 


1) Cooperation in Spectrum Sensing: Modeling and perfor¬ 
mance optimization of cooperative sensing schemes has long 
been of great interest. Several cooperation methods have been 
introduced and analyzed in the literature, such as the scheme 
based on voting rules |26) as a simple suboptimal solution, or 
the so-called OR logic operation )27) and AND logic operation 
[28). The likelihood-ratio test (LRT) )29) , (30) is known as 
the optimal fusion method when the distributed nodes report 
their sensing outcomes to the FC through nonideal analogue 
communication links The effect of reporting channel 

impairments on the overall sensing performance has been 
investigated in [5), pT)-|33~). Particularly, the authors in [311 
compare the performances of the hard-decision- and soft- 
decision-based fusion methods and illustrate that, in general, 
the soft decision significantly outperforms the hard decision 
when nonideal reporting channels are considered. Although 
LRT is commonly considered as the best soft-decision-based 
cooperation, finding the optimal LRT thresholds for individ¬ 
ual nodes and for the FC incurs prohibitive computational 
complexity ®, pT) . Consequently, optimal linear combining 
©-0 has been suggested as a very good alternative which 
provides nearly-optimal results but at affordable computational 
cost. As a fast alternative linear combining approach with a 
very good performance, the so-called deflection criterion [ |29) , 
[34), a.k.a., deflection coefficient (DC) |35| or its modified 
version, modified deflection coefficient (MDC) 0, ®, Q5), 
is commonly used in the literature to design linear 
fusion schemes in cooperative sensing scenarios. Moreover, 
the definition of DC is extended in 0 and Hz) to enable joint 
optimization of reporting and fusion processes in cooperative 
spectrum sensing. 


2) Censoring in Cooperative Sensing: Censoring-based dis¬ 
tributed detection is investigated in m where the system 
model is provided both in the Bayesian and Neyman-Pearson 
settings, a so-called censored likelihood ratio test is defined, 
and censoring regions are determined for a certain overall de¬ 
tection performance. In [38), performance of OR-based coop¬ 
erative sensing with censoring is optimized subject to an upper 
bound on the cost of local sensing and reporting processes. 
The results provided in [9;| show that the communication 
resources required for reporting is significantly reduced when 
a censoring-based cooperation is established between the sens¬ 
ing nodes. Efficient censoring-based cooperative sensing under 
bandwidth constraint is discussed in [39] where performance 
of hard-decision based spectmm sensing is studied with both 
perfect and imperfect reporting channels. In |]Tl), analytical 
expressions for the sensing parameters are derived according to 
a Neyman-Pearson setup for censoring with both soft and hard 
fusion schemes, but no constraint on the energy consumption 
is taken into account. In © censoring in a wireless sensor 
network is considered where all sensors experience the same 
SNR and have the same performance while TDMA is assumed 
for reporting the local decisions to the FC through ideal, i.e., 
error-free channels. Based on this model, the throughput is 
maximized for a k-out-of-N fusion rule under a constraint on 
energy consumption. The censoring effect in hard-decision- 


based cooperative sensing is investigated in © with two 
different sensing methods. In the first method, the number of 
samples used by each sensor is fixed whereas in the second 
approach a so-called sequential sensing technique is used 
where the sensors sequentially collect observation samples 
until they reach a decision about the presence or absence of the 
PU. The latter approach is shown to be more energy efficient 
than the former one for a given performance. A combined 
sleeping and censoring method is proposed in © which 
serves as an energy efficient technique for cognitive sensor 
networks with hard decision fusion. The objective in © is to 
minimize the energy consumed in distributed sensing with OR 
fusion mle subject to constraints on the detection performance, 
by optimally choosing the sleeping and censoring design 
parameters. The constraint on the detection performance is 
given by a minimum target probability of detection and a 
maximum permissible probability of false alarm while the 
local SNRs at the sensing nodes are assumed to be the same for 
analytical tractability. The combined censoring and sleeping 
approach is further investigated in © where sensing param¬ 
eters are derived by minimizing the maximum average energy 
consumption per sensor subject to a lower-bound on the overall 
probability of detection and an upper-bound on the overall 
probability of false alarm for both OR and AND fusion rules. 
Compared with [ 151 which intends to minimizes the network 
energy consumption, in © minimizing the maximum average 
energy consumption per sensor is considered. 

3) Sensor Selection/Scheduling Algorithms: In ©, © 
sensor selection methods are developed to increase the energy 
efficiency of cooperative sensing in CRNs with hard decision 
fusion. They concentrate on finding one subset that has the 
minimum energy consumption. A sensor selection problem for 
cognitive wireless sensor networks with hard decision fusion 
rule and error-free reporting channels is considered in (19 | 
where the focus is on maximizing the network lifetime. In 1201 
an energy-based sensor selection algorithm is considered in a 
cognitive wireless sensor network to provide approximately 
the same lifetimes for sensors via appropriate design of 
cooperative sensing with hard decision fusion. It is assumed 
in pO) that all of the sensors experience the same SNR and 
operate with the same sensing thresholds. An energy efficient 
sensor scheduling algorithm is proposed in pT[ to optimally 
schedule the activities of sensors to provide the required 
sensing performance and increase the overall secondary system 
throughput. This sensor scheduling method is based on OR 
fusion rule and assumes error-free reporting channels. Sensor 
scheduling schemes in |22)-)25) aim at extending the network 
lifetime while maintaining full coverage in a certain area. 


B. Contribution 

In this paper, we extend the works in [6)-[® and propose a 
novel energy-efficient cooperative sensing scheme for CRNs. 
This new design is achieved by adding two new mechanisms 
to the commonly-used cooperative sensing structure. The first 
mechanism is realized as a set of controllers mounted on the 
sensing nodes to make random energy-saving interruptions in 
their cooperation with the FC. The second mechanism is a 
compensation process at the FC. This compensation, which is 
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realized as a linear estimator, aims at recovering the local test 
summaries out of degradations caused jointly by the interrup¬ 
tions and reporting channel contaminations. The estimation of 
the local test summaries is realized in the proposed system 
by using the spatio-temporal cross-correlations of the sensor 
outcomes as well as auto-covariance functions characterizing 
the behavior of the reporting channels. 

We model and thoroughly analyze the proposed system to 
derive the global test summary statistics in terms of probability 
distributions based on which the controllers work. By using 
the statistics obtained, we formulate the proposed system 
performance optimization according to two commonly-used 
criteria, namely the Neyman-Pearson criterion (NPC) and the 
deflection criterion. Specifically, first, we analyze the global 
test summary statistics and characterize the system perfor¬ 
mance by providing closed-form relations for the probability 
of false alarm and the probability of detection. We then derive 
the detection threshold for a fixed false alarm probability and 
formulate the system performance optimization based on NPC 
and subject to a constraint on the energy consumed at the local¬ 
sensing and reporting phases. In consequence, we formulate 
a significant tradeoff concerning the overall detection quality 
along with joint energy consumption of the local sensing and 
reporting phases. We discuss that the obtained optimization 
is inherently an stochastic program with prohibitive compu¬ 
tational complexity and then, we consider the optimization 
as a decision making process with uncertainty. By using this 
decision making structure, we show the link between the NPC 
and DC for the proposed cooperative sensing scheme. 

Regarding the DC, we formulate the performance optimiza¬ 
tion as a convex-over-convex fractional programming problem. 
This fractional program, which is solved by a branch-and- 
bound (BnB) algorithm, gives optimal distributions of random 
interruptions for maximizing the overall detection perfor¬ 
mance, subject to a constraint on the energy consumed at the 
local-sensing and reporting phases. Through delicate algebraic 
manipulations, we convert the fractional program into a set 
of nonconvex quadratically-constrained quadratic programs 
(QCQP) whose solution can be obtained in polynomial time 
by using standard techniques in nonconvex quadratic program¬ 
ming. Therefore, we construct an effective computationally- 
affordable method for solving the proposed optimization even 
when dealing with large-scale networks. 

C. Discussion and Comparison 

In the proposed method the local sensing and reporting 
processes in a CR node are occasionally suspended if the con¬ 
tribution of that node to the overall sensing process does not 
significantly improve the system performance. As explained in 
the following, the result of such an interruption is more than 
just saving the energy of the sensing nodes. 

Spectrum sensing in CRNs is often concerned with protect¬ 
ing several, and in many cases, heterogeneous PUs. Moreover, 
the sensing process is often performed on multiple frequency 
bands. When a CR node is not able to significantly contribute 
to the overall sensing process, this generally means that one 
or maybe more PUs operating on a certain group of frequency 
bands are hidden from this particular sensing node. However, 


this CR node may still be able to detect signals from other PUs 
operating possibly on a different geographical location and/or 
on a different frequency band. Therefore, instead of perform¬ 
ing spectrum sensing continuously and reporting unreliable 
results which will most-likely be discarded at the FC, it is more 
efficient for this CR node not to waste its resources in sensing 
the signals from the hidden PUs and instead, to allocate its 
signal processing and reporting resources to discover signals 
from other PUs which are potentially not hidden. 

Therefore, by using the proposed cooperation scheme, the 
sensing nodes can spend their available resources more effec¬ 
tively and obtain a better cooperation performance based on 
the available resources. Note also that, this efficiency cannot 
be achieved by the optimal linear cooperation in which the 
sensing nodes continuously perform their spectrum sensing 
and it is up to the FC to decide about their contribution by 
assigning weights to the received reports. 

The censoring-based methods are not as effective as the 
proposed random interruptions either. This is due to the fact 
that, in cooperative sensing with censoring, a sensing node 
continuously performs the local sensing process regardless 
of the quality of the PU signal at its receiver, and only 
the reporting process is avoided when the sensing node is 
not able to reliably detect the PU signal. Therefore, in a 
censoring-based cooperation the CR nodes with unreliable 
sensing outcomes are not able to save their local sensing 
resources or use them to sense other PU signals or other 
spectrum bands. 

In sensor selection or sensor scheduling algorithms, energy 
efficiency is achieved by giving binary send/no-send instruc¬ 
tions to the sensing nodes. Despite the variety of objective 
functions used in their performance optimization, these al¬ 
gorithms can be considered as a special case compared to 
the proposed interruption-based cooperation. The reason is 
that even though the efficiency is realized in the proposed 
method based on binary sleep/wakeup instructions, the overall 
contribution of each node is not affected in a binary manner. 
As a simplified explanation (see Section |Il-B| and Section [Hi] 
for details), since we are dealing with the statistics of the 
received reports which go through a compensator and a soft 
fusion process, the contribution of each node (let’s say node 
i), is affected (i.e., multiplied) by a continuous factor (i.e., 
the Bernoulli probability pi). Therefore, the reporting power 
and local computation cost are scaled down by a continuous 
factor, compared to a simple send/no-send structure which 
turns the sensing node either on or off. Nevertheless, the 
proposed method still enjoys the benefits of an on-off power 
saving scheme due to its unique sleep/wakeup structure. 

To the best of our knowledge, the present work is the first 
one to consider energy efficiency in a CRN by optimizing the 
detection performance of cooperative spectrum sensing with 
soft fusion rule at the FC, nonideal reporting channels, and 
different SNR levels at the sensing nodes which do not censor 
their local sensing outcomes. 

Organization 

The rest of this paper is organized as follows. In Section 
[II] general modeling assumptions and details of the proposed 




4 


system are introduced and its behavior is analyzed. In Section 
PH the system performance optimization is formulated in a 
Neyman-Pearson setting and the challenges in the optimization 
are discussed. In Section [TV] a scenario-tree-based stochastic 
programming method is developed for the problem. The 
stochastic programming approach is then linked to the DC in 
Section [V] where a low-complexity effective solution for the 
problem based on standard semidefinite programming methods 
is constructed. The effectiveness of the proposed cooperative 
spectrum sensing is then demonstrated through simulation re¬ 
sults in Section VI followed by concluding remarks provided 
in Section IVTlI 


II. System Model and Analysis 
In this section, we first introduce our notation along with a 
set of definitions which facilitate mathematical representation 
of the proposed analysis. Then, the system model is intro¬ 
duced in detail and the parameters needed for performance 
optimization are derived accordingly. 




= E 


= R, 


( X -Mx|«j(y- aviwJ Uh] 


H 


:y| 'H-h l 1 x\'H h ^y\'H h 


( 6 ) 


where and ^ y \-n h denote mean of x and mean of 

y, conditioned on Hh, respectively. The main diagonal of 
the autocorrelation matrix R x |'H h is denoted CK x |-^ h which 
is defined as 


(IH x|h J. ,^^(R x | Hh ). . (7) 

where 6,, :) is Kronecker’s delta function which equals to 1 if 
i = j and 0 otherwise. 

Non-conditional statistics are related to their conditional 
counterparts according to the following definitions 

C x = Pr{7f 0 }C x |„ o + Pt{H r}C x|Wl (8) 


C xy = Pr{7f 0 }C xy |„ 0 + Pr{7fi}C xy | Wl (9) 


A. Notation 

Matrices and column-vectors are denoted in boldface by 
uppercase and lowercase letters, respectively. The notation ||x|| 
denotes the Euclidean norm of the vector x. The notation 
diag(x) represents a diagonal matrix whose main diagonal 
is x. Matrix inequality is represented by A, i.e., A A B 
means that A B is positive semidefinite and for a vector 
it denotes the element-wise inequality with A representing 
the strict inequality. The identity matrix is denoted by I and 
0 n denotes the n x n null matrix while 1 and 0 denote all- 
ones and all-zeros column vectors, respectively, e, denotes a 
column vector where all elements are zero except for the ith 
element which is one. Tr(A) refers to the trace of A, i.e., 
sum of elements in the main diagonal of A. The vectorization 
of matrix A is denoted by vec (A) which represents a vector 
obtained by stacking the columns of A on top of one another. 
AoB represents the Hadamard product of matrices A and B. 

In order to account for temporal and spatial representations 
of signals, we use the following notation. For vector x in the 
following form 

x(m) = [xi(m), ...,x K {m)] T (1) 

we consider these three notations 

X L (m) = [x(m),x(m - 1), ...,x(m - L)] T (2) 
Xi(rn) = vec (X L (m)) (3) 

X L (m) = diag (x L (m)) (4) 


Note that we represent the time index by m. 

For the two hypotheses considered in this paper, i.e., the 
null hypothesis Ho corresponding to the absence of the PU 
signal and H\ representing the presence of the PU signal, the 
following notations are used to represent conditional second- 
order statistics of signals (for h = 1,2) 


"'x| Hh 


= E 


= R» 


( x ~ /^x|«J( x - MxlwJ ^h\ 


H 


I 'Hh ( J ‘x\'H. h f J ‘x\'Hh 


(5) 


Note that other non-conditional statistics in this paper are 
related to their conditional counterparts in a similar fashion. 

As another definition regarding the second-order statistics, 
we have the following K{L + 1) x K matrix 


C xty = E 


x L(m)y\ (to), ..., xl (m)y' K (m) 



( 10 ) 


where // x . = E[xi] and f denotes complex conjugation. 

Scalar conditional autocorrelation and autocovariance func¬ 
tions of x(m) for k = 1,2,..., A, n = 1,2, l = 

0,1,..., L, r = 0,1,..., L, and h = 0,1 are denoted by 

w x | H h (k,n;l,r) = E [x k (m - l)xl(m - r)[H h ] (11) 


Cx|H h (k, n; l, r) = w x | Hh (k, n; /, r) 

- E [x k (m - l)\Hh ] E [a;* (to - r)\H h \ 

( 12 ) 

Note that with this definition for the autocovariance function 
we have (for i = 1,..., K and j = 1,..., K ) 

(CxiwJ^j ; = c x \ Hh (iJ-,0,0) (13) 

However, regarding the covariance matrix of x.l, which is 
denoted by C Xjj | Wfe , we have (for i = 1 ,...,K(L + 1) and 
j = 1,..., K(L + 1)) 


(C XL \n h )ij = Cx\u h (k,n;l,r) (14) 


where k = 


L+l 


n = 


L+l 


, l = (i + L)mod(L + 1), 


and r = (j + L)mod(L + 1). [•] and mod denote the ceiling 
function and modulo operation, respectively. 

Now we are ready to model and analyze the proposed 
system behavior. 
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Fusion Center 


Fig. 1. The proposed system architecture. 



Fig. 2. Controlling the local sensing and reporting processes at the sensing 
nodes based on the predetermined sequence of binary random numbers which 
are generated and shared with CR nodes when initializing the cooperation. 

B. System Model 

A CRN with K sensing nodes is considered. These nodes 
cooperatively sense the radio spectrum to find temporal and/or 
spatial vacant bands for their data communication. Fig. |T] 
shows the major elements of the proposed design. Each CR 
node is equipped with a built-in spectrum sensor which enables 
it to detect the PU signal through inspecting its own listening 
channel. Listening channels are referred to the channels be¬ 
tween the PU and the sensing nodes. The sensing nodes have 
access to a dedicated but nonideal reporting channel to send 
their individual sensing outcomes to the FC. 

In our adopted model, the mth sample of the received PU 
signal at the ith CR node is represented as 


assumed to be in the form of energy detection which can be 
represented as 

IV—1 

Ui(m) = E \xi(Nm — k)\ 2 (16) 

k—0 

The resulting sensing outcome is then transmitted to the FC 
through the reporting channel. We represent this interrupted 
cooperation as 

yi{m) = 9i(m)ui(m ) + Zi(m) (17) 

where yi(m) denotes the ith received sensing outcome at 
the FC, 9i{m) £ {0,1} denotes the random number which 
controls the zth sensing node, and Zi(m) ~ C7V(0, <r^) is 
the reporting channel contamination modeled as AWGN. It is 
worth noting that, for simplicity, we have not explicitly shown 
the controllers in Fig. |T| However, © clearly demonstrates 
that when the generated random number 0, is zero, nothing 
is transmitted to the FC. This process can be implemented in 
practice by using an energy-saving controller which suspends 
the local sensing and reporting processes when the correspond¬ 
ing random number is zero, see Fig. [2] 

In the matrix form we can represent ( flTj i as 

y(m) = diag(0(m))u(m) + z(m) (18) 

where y(m) = [yi(m), y K {rn)] T , 0(m^ = 

[0i {m),...,6 K {m)] T , u (m) = [ui (m),...,«jr(m)] T , and 

z( to) = ..., ZK{tn)\ ■ In the proposed system the 

random numbers 9 t , i = 1, K, are generated independently 
over time and space, i.e., 9 k {m — l) and 9 n (m — r) are 
independent if k ^ n or / ^ r. Moreover, 9i(m) follows 
the Bernoulli distribution with Pi(m) = Pr {9i{m) = 1}. 
Therefore, the autocorrelation and autocovariance functions 
of 6 can be represented as 

uig{k,n;l,r ) = p k (m - l)p n (m-r ) 

+ S k , n Si,rPk{m - /)[ 1 — Pk(m - l )] (19) 


f Xi(m) = Ho 

\ Xi(m) = his(m) + Hi 

where s(m) denotes the signal transmitted by the PU and 
Xi{m) is the received signal by the zth SU. hi is the listening 
channel block fading coefficient, which is assumed to be 
constant during the detection interval. Listening channel gains 
are assumed to be independent circularly-symmetric Gaussian 
random variables. Vi(in) denotes the circularly-symmetric 
zero-mean additive white Gaussian noise (AWGN) at the CR 
sensor receiver, i.e., zzj(ro) ~ CAT(0, a ^.). s(m) and 
are assumed to be independent of each other. 

In this cooperative sensing scheme, the behavior of each 
node is controlled by a predetermined sequence of binary ran¬ 
dom numbers generated at the FC. Whenever the cooperative 
sensing is performed, only the nodes whose corresponding 
random number is one are allowed to contribute to the overall 
sensing process. Specifically, if the zth CR is among the 
cooperating nodes, i £ {1, 2, K }, it performs spectrum 
sensing by using its built-in sensor to derive a local test statistic 
Ui. Without loss of generality, the local sensing process is 


cg(k, n; l, r ) = u g {k, n; l, r ) - p k {m - l)p n {m - r ) 

= Sk, n Si,rPk{m — Z)[l — Pk{m - /)] (20) 

Moreover, the reporting channel is characterized by 

(k, n; l, r ) = <j 2 z 5 k ,n8i, r ( 21 ) 

Note that 9, u, and z are assumed to be independent of 
each other. We also need p( to) = [pi(m), ...,pif(?7i)] T and 
b(m) = in our formulations. b(m) de¬ 

notes a realization of the random vector 6{m). For simplicity, 
we will drop the time index m when representing vectors and 
matrices. 

The reported sensing outcomes feed a linear estimator, 
depicted in Fig. [T]as compensator, which uses spatial and tem¬ 
poral correlations of the received test summaries to estimate 
their actual (i.e., non-contaminated and non-interrupted) val¬ 
ues. Specifically, the observation vector used by the estimator 
is y l which, based on the notation provided in Section |II-A| 
can be expressed as 


y L = ©lUl + z; L 


( 22 ) 
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Note that @a is a K(L + 1) x K(L + 1) diagonal matrix. The 
estimation process at the FC is represented as 


u = £, T y L 


(23) 


where the weight vector £ and the constant e are obtained as 
the minimizer of the mean-squared error (MMSE) 


follow the Gaussian distribution. Consequently, the estimated 
local test summaries conditioned on a given realization of Ol 
are normal random variables, i.e., their conditional probability 
density function (pdf) can be represented as (for h = 0,1) 

f (u\n h ,e L = b L ) = 


J(S,e)=E 

l|b-u|| 2 

= E 

£ T y l + e - u 

2 

■ (24) 

1 



- 



x exp 

~ 9 


(2n) K \C u]iHM \ 


Accordingly, the MMSE weight vector can be obtained as 

r = C-}C yLU . ( 25 ) 


(*■ 




i| {-H, 




t J, u\{'Hh,b L } 

(31) 


where 


Moreover, since u has nonzero mean in general, the optimum 
value of the bias term e is given by 

e* = E[u] - C T E\y L ]. (26) 


The use of linear estimation—instead of e.g., maximum 
a posteriori (MAP) estimation—in the proposed design is 
supported by several factors. First, to implement a linear 
estimator only first- and second-order statistics of the local 
sensing outcomes are required. Second, the analysis based 
on linear estimation leads to tractable formulations which 
facilitate performance optimization. And third, as it is shown 
in Section [Vi] the proposed structure achieves the optimal 
performance of existing cooperative sensing schemes but with 
a significant level of energy efficiency. 

From ( |25) l and ( [26] ) we see that in order to establish the 
optimal linear estimation, we need C yj U and C y; . We derive 
these covariances in terms of the statistics of the local sensing 
outcomes and random interruptions as follows. First note that 
Cy LU is obtained as 


Cy tU - PiC UtU . 


(27) 


Note that = E\&^]. Then, by using (19 1 —(21 1 the elements 
of C yx i^ h can be obtained as in (281 at the top of the 

following page. Recall from Section II-A that k = 


L+1 


jq-j- , l = (* + L)mod(L + l), and r = (j + L)mod(L + 


n = 

1). Hence, we have shown that 


Cy L \n h =PL (i-P l) 
and consequently. 


^u L | U h + Pi C u L |HfcPi + 

(29) 


C yL = P L (i - P L ) 5 R Ul + PlC Ul P l + <7*1. (30) 

It is worth noting that, C U(U , C UL , and 91 u , are obtained 
in terms of the second-order statistics of the channel gains 
between the PU and SUs as well as noise variances experi¬ 
enced by the sensors. The details are provided in Appendix I. 
Availability of these channel statistics is a common assumption 
in designing cooperative sensing schemes for CRNs §-§• 
Since most PUs of interest such as current TV stations and 
OFDM systems operate based on certain standards, transmit at 
fixed power levels, and/or periodically transmit pilot signals, 
even exact channel gains can be obtained in many practical 
cases, see e.g., Ill- 

According to the central limit theorem, if the number of 
samples N in ([T6]l is large enough, the local test summaries 


Vb\{H h ,b L } - E [u| Kh, 0 L = b L ] 


(32) 


= C E @lUl + z L \T-Lh, Ol = bz 

= C T B L E[u L \H h ]+e* 


and Cu| {Wh . bl/ } denotes the conditional covariance matrix of 
u. Note that, conditioned on 0/ = b^, the observation vector 
is yL = BlUl + zl which leads to 


C yL\{H h .b L } = BLC Ut | Wh B z 


(33) 


Therefore, by using (23 i we have 


c» mh , bL} =e T c yLl{Hh , bL} e 

= e T (B L c Ui |„ h B L + c ZL )r 

= e T (B L C ULlHh B L +a 2 z l) r 


(34) 


Now, the total probability theorem gives the distribution of the 
estimated local test summaries, i.e.. 


2 n — l 

/ (u| Uh) = Y, /a ("I U h , e L = h L ) Pr {0 L = b L } (35) 

b_L=0 


where n = K(L + 1) and 


Pr {0 L = b^} 

L K 

=n n p^ m - i) bk(m ~ i) [i -Pk{m~ i 

;=0 fc=l 

(36) 

Therefore, for a given we can establish the optimal fusion 
rule at the FC as the well-known FRT, i.e., 




/Ml> 

f(u\n 0 )</ 


(37) 


where the threshold ta is determined based on minimizing the 
average Bayesian cost of detection. 

It is worth pointing out that, in existing cooperative sens¬ 
ing schemes J6l-|8l where there is no interruption nor is 
any compensatior0 the optimal performance of the FRT 
detection is closely achieved by a linear detector where the 
weighting vector is derived in terms of the listening-channel 
and reporting-channel statistics. More specifically, a so-called 
global test summary, which is used to decide the presence or 


'Note that in this case u = u + z. 
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= E l( 6 k{m - l)u k (m - l) + z k (m - l )) (0„(m - r)u n (m - r) + z n {m - r)) \H h \ 

- E [( 8 k (m - l)u k (m - l)) \ T-L h \ E [( 9 n (m - r)u n {m - r)) \H h \ 

= C 0 (k,n\l,r)u} u \- Hh (k,n;l,r) +p k (m - l)p n (m - r)c u \ Hh (k,n-,l,r) + w z (k,n;l,r) 

= S k , n 5i, r Pk(m - /)(1 — Pk(rn - l))uj u \ nh (k,k;l,l) + p k (m - l)p n {m - r)c u \ nh (k,n;l,r) + a 2 ,5 kjn 8i tr 
= 8 i,jPk{m - 0(1 -Pk(rn - l))co u \ nh (k,k;l,l) + p k (m - l)p n (m - r)c u |^ h (fc, n; l, r) + cr^- 
= P L (i- P L ) 91 UL | Wh +P L C UL | Wh P L + all 

L \ / J i,j 


(28) 


absence of the PU signal, is constructed as a weighted sum 
of the received reported sensing outcomes. In the proposed 
system, the linear combining is represented as 

Ui 

S'(u) = w 7 u ^ rs (38) 

n 0 

where S(u) is the global test summary. 

The main idea of linear combining is that the combining 
weight for the signal from a particular user represents its 
contribution to the global decision. For example, if a CR 
generates a reliable (i.e., high-SNR) signal that may lead to 
correct detection on its own, it should be assigned a larger 
weighting coefficient. For those SUs experiencing deep fading 
or shadowing, their weights are decreased in order to reduce 
their negative contribution to the decision fusion. 

Although linear combining is known to be an effective and 
computationally-affordable technique in designing cooperative 
sensing schemes, it is not energy-efficient. It is worth noting 
that, regardless of the channel conditions experienced by dif¬ 
ferent sensing nodes, they all provide their sensing outcomes at 
the same cost. This cost manifests itself in the form of energy 
consumed at the local sensing and reporting phases. When 
the contribution of a particular node to the global decision 
is suppressed by assigning a small weight to it, the energy 
consumed by that node for the local sensing and reporting 
processes is relatively wasted. 

In the proposed design, instead of suppressing the contri¬ 
bution of nodes working under deep fading or shadowing, 
they are occasionally ordered not to cooperate. In this struc¬ 
ture, the more reliable a sensor is, the more likely it is to 
contribute to the overall sensing process. Consequently, the 
required discrimination between the reliable and unreliable 
sensors is achieved while no energy is wasted. Moreover, it is 
possible to optimize this process considering a constraint on 
the energy consumed at the local sensing and reporting phases. 
More specifically, we are interested in the optimal distribution 
associated with the random number generators, i.e., optimal 
value of pi, assuming a general linear fusion rule at the FC. 

Remark 1 : The random numbers are not required to be sent 
to the CR nodes during the sensing process. Hence, there is 
no intervention required from the FC while the local sensing 
is performed by the CR nodes. In fact, communication of 
the random numbers takes place before starting the sensing 
process. Based on the listening-channel and reporting-channel 
statistics, the distribution of the binary random numbers are 


first calculated at the FC and then, each CR node can 
either receive a random bit string generated according to 
the corresponding Bernoulli distribution, or receive only two 
numbers to feed a simple binary random number generator. 
One of the numbers indicates the corresponding Bernoulli 
distribution and the other one serves as a seed for the local 
random number generator. Once initialized, the cooperative 
sensing is performed in a distributed manner without any 
intervention from the FC. We know that in practical CRN 
architectures, such as the one assumed in the IEEE 802.22 
standard the FC regularly communicates with the CR 
nodes to coordinate their communication-related activities as 
well as their local sensing and reporting processes. Therefore, 
the proposed cooperation structure fits easily to the commonly- 
used CRN configurations. 

III. Problem Formulation 

In binary hypothesis testing for signal detection, the detector 
performance is commonly characterized by two probabilities, 
namely, the probability of false alarm Pf and the probability of 
detection P&. The system performance optimization based on 
NPC requires maximizing the probability of detection subject 
to a constant false alarm probability, i.e., 

max Pd s.t. Pf = a (39) 

where a is referred to as the significance level of the detector. 

Detection and false alarm probabilities depend on the global 
test summary statistics. Since and z k follow normal 
distributions, we can conclude from ( |22| ), ( |23j ), and ( |38| that, 
conditioned on 0 k = b /, our test summary ,S'(u) is a normal 
random variable with the following mean and variance (for 
ft = 0 , 1 ) 

E[S{u)\'H h ,6 L = b L \ = w T /r u|{ « h , bt} (40) 

Var [S(u)\n h , 0 L = b L ] = w T C u | {WhjbL} w (41) 

Consequently, the detection and false alarm probabilities, 
conditioned on 0 k =h k , can be expressed as 

P d \ hL =Pr{S>T S \n 1 ,e L = b L } (42) 

TS- W T ^ |{ ^ 1|b ^ } \ 

y / w T c aKWl , bi} w j 







P i \ hL =Pr{S>rs\no,e L = b L } (43) 

= Q ( Ts - wT ^\{n 0 ,b L } 

V Y / w T C ii | {Woibi} w 

where P f | bi and / d | b; denote the false alarm and detection 
probabilities conditioned on 6l = t>L, respectively. The total 
probability theorem then gives the overall detection and false 
alarm probabilities based on their conditional counterparts, i.e., 

2"-l 

Pd= 51 P d| bi Pr{0L = M (44) 

h L —0 
2 n — 1 

Pf = E P f|b,Pr{0L =b L } (45) 

bi ,=0 


Given the target false alarm probability Pf = a, the 
detection probability can be obtained by removing rg from 
© and ©■ More specifically, solving Pf\b L = a for T s 
gives 

TS = Q _ 1 (a)y / w T C b | {Woibi} w + w^VuKWo.b,,} ( 46 ) 

According to ( |46| , having a constant false alarm probability 
requires a detection threshold which depends on the random 
number generator outcomes. In other words, unlike in 
®> the detection threshold for a given false alarm probability 
in the proposed structure is not fixed and in fact, can be 
viewed as a random variable determined as a function of 0 [ . 
The proposed threshold adaptation is similar to the adaptation 
process in radar receivers with constant false alarm proba¬ 
bility, a.k.a., constant false alarm rate (CFAR) systems. The 
difference here is that the threshold depends not only on the 
received noise level, but also on the set of nodes participating 
in the cooperative sensing. 

The detection threshold ts, when plugged into ( |42| , leads 
to the following conditional detection probability 


p(“) 

Gilbi 


= Q 


Q (a)V wTc a|{w 0 ,b i }W-a^ >L 

y / w T C 0l | {WljbL} w 


(47) 


( a ) 

where P^| bt denotes the conditional detection probability 
when the system false alarm probability is set to a, and 


a b L 


A_ 

~ A t u|{'Hi,b L } ^ulJHo.bi} 

= r T B l (E [u l |^i]-S[u l |H 0 ]) 


= C 


T 

U L U 


P L 



XB L (E [u L |Hi]-£7[u L |fto]) 


PlC U l P l + ol I 

(48) 


Moreover, by using <[25]>, <[29|), and ([34]) we have 


i|{«h,bi} 


C^ u Pl 




^ Ul +PlC Ul P, 


ail 


i -i 


x 

X 


(B L C Ut| « h B L+( 7 z 2 l) 

Pl (i - Pl) ^H Ul + PlC U l P l + o\l 


Pr.C, 


( 49 ) 


In order to clearly see the role of in the proposed system 
performance, we use the fact that the elements of P / lie 
between zero and one. Accordingly, we apply approximations 
on ( |48] > and ( |49[ i based on first-order Taylor series expansion, 
i.e., (I + A) - « I — A when A fe —) 0„ as k grows. As the 
result, we obtain 

a bi « CI lU P l B l (E [u l |Hi] - E [u L |H 0 ]) (50) 

c a|{« h ,bi} ~ C u tU P l (B L C Ul/ | Wh B L + a'llj P l C U lU 

(51) 


Through some further algebraic manipulations, the condi¬ 
tional detection probability is reformulated as 


p(«) 

a|b L 


(Pi) ~ Q 


where. 


Q 1 (ot)Jpl^ {Ho ,h L }PL-q L l L PL 


,bi}PL 


(52) 


q bi = C utu w o B l {E [u l |Hi] - E [u l \H 0 ]) (53) 
V {HM 4 (B L C ULlHh B L +l) o (C UiU ww T Cj iU ) 

(54) 


For simplicity, we have assumed that cr 2 = 1. It can be easily 
shown that bi j is positive semidefinite and the elements 

of q b( are all nonnegative. 

From ( [52] ) the overall detection probability for the target 
false alarm probability of a can be obtained as 

Pt\pL) = E P = M (55) 

b i= 0 

(a) 

where P d denotes the overall detection probability when the 
system exhibits a constant false alarm probability of a. 

Based on the analysis provided, now we can represent our 
proposed detector performance optimization as 

Pl = argmaxP d (Q) (p L ) (Pl) 

PL 

s.t. 0 A p L A 1 , l 7 • pl < (1 - r))n 


where rj denotes the power-efficiency constraint we consider 
in designing the proposed interrupted reporting scheme. Note 
that for 0 < 77 < 1 some elements of pl have to be less 
than 1 , which means we are forcing some sensing nodes to 
occasionally go to the sleep mode (i.e., avoid cooperation). 
Moreover, it is clear that in general, the higher 77 we chose, 
the higher is the number of nodes which are forced to sleep, or 
the more likely any node is to go to the sleep mode. In fact, 77 
is the parameter by which we control the energy consumption 
of the proposed interrupted reporting scheme. 

Remark 2: Compared to the cooperation method in |hj [ 8 |. 
the proposed approach consumes less energy. Specifically, in 
0 -® before starting the sensing process, a set of weights are 
calculated to discriminate between the sensing nodes whereas 
in the proposed method a set of Bernoulli probabilities are 
calculated with the same objective. Since we have not changed 
any of the major phases in the general cooperation process 
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( = 0 


t = \ 


t = 2 



Fig. 3. Scenario tree for the proposed T-stage stochastic optimization. 


considered in those works, it can be concluded that, con¬ 
sidering the entire sensing process, the proposed cooperation 
scheme is more efficient. 

Our optimization in is a challenging problem since, 
firstly, as shown by ( [55] ), the objective function is constructed 
as a summation over 2” nonlinear functions of p l, and sec¬ 
ondly, the summand functions are not concave. Solving such 
an optimization which needs exponentially-growing number 
of operations to be built would be a serious challenge even 
if all the summand functions were concave. In the following 
section, we tackle these issues by developing a stochastic 
optimization approach for the problem through considering 
€D as a decision making process with uncertainty. This 
approach will lead us to construct an effective low-complexity 
solution for the performance optimization of the proposed 
system in Section [V] 

IV. Stochastic Optimization 

In this section, we develop a multistage stochastic program¬ 
ming approach based on a so-called scenario tree structure 
ED- We first use this scenario tree to represent an alternative 
formulation of the optimization problem. This new formula¬ 
tion enables us to derive an alternative problem formulation 
whose objective function is built more conveniently. Next, we 
exploit the scenario tree to construct a two-stage stochastic 
optimization scheme which can serve as a fast procedure 
in the proposed system to obtain an approximate solution. 
The decision making structure developed in this section along 
with the two-stage stochastic optimization are used in Section 
|V] to develop, based on standard semidefinite programming 
techniques, a polynomial-complexity optimization approach 
which gives very good results. 

First note that the objective function in €D can be viewed 
as the expectation of the performance metric (i.e., the prob¬ 
ability of detection) over the probability distribution of 6 l. 


i.e.. 


P^Pl) 


e 0 l 



(56) 


From the stochastic programming point of view, 0 / is 
the source of uncertainty in our decision making process. 
Hence, successive realizations of 0 r form the branches of the 
scenario-tree considered, see Fig [3] 

The proposed scenario tree is constructed as follows. A root 
node is associated to the first decision stage and to the initial 
absence of observations. The first decision stage is designated 
as t = 0. The nodes associated with decision stage t, t = 
1,..., T, are connected as children nodes to the nodes in stage 
t — 1. The observation at the decision stage t is denoted by 
ip t and defined as 


ip t = h L (m + (t- !)(£ + 1)) (57) 


Each node at stage t corresponds to one possible outcome 
for ip t given the observation corresponding to its parent, 
i.e., tp t _ 1 . The observations associated to the nodes on the 
unique path from the root to a leaf define together a particular 
scenario. Specifically, the fcth scenario is denoted by i]) {k> = 
[tp ^,..., which is associated with the fcth branch of the 

scenario tree while ip t represents the observation on the fcth 
branch at stage t. 

(k) 

At stage t of the fcth branch and after observing xp t , 

(k) 

a decision is made, denoted by Q ’, which determines the 
Bernoulli distribution of the random variables for the next 
stage, i.e., for t = 0,1, we set 

p L (m + t(L + l)) := Ct k) (58) 


Recall that pi(m) = Pr { 6 ^( 171 ) = bi(m)}. Moreover, we 
know that each realization of 0 L determines the state of the 
sensing nodes for L + 1 consecutive time slots. Hence, the 
indexing distance between two consecutive decision stages 
should be (L + 1) as indicated in ( [57] ) and ( |58] >. Note also 
that, the initial decision is made with no observations 
available. 

We collect all the decisions associated with the fcth scenario 
in = [Co*^ Ci fe \ Ct )]• B ase d on these decisions, the 
occurrence probability of the fcth scenario, denoted by n lk> , 
is derived as 


7 t (ft) = Pr{i/> (fe) } 

= Pr{t/>! fc) } Pr • ■ -Pr|^r Vt-d‘ " >^ 1 *°} 

= <P V (V4 fc \ C^) (Vr\ C { t-i) (59) 

where ip (tp[ k \ Ct-i) denotes the probability of observing 

(k) 

ip t at stage t. This probability is calculated based on decided 
Bernoulli distributions as follows 



=nc-iw 


np/’t (* 


1 - cK« 


-,1 


(60) 


where index i refers to the vth element of the vectors. 

Considering the fcth scenario, based on the value of the 

(k) 

random vector realized at each stage, i.e., ’, and Bernoulli 



















10 


(k) 

distributions assigned, i.e., Q , the probability of detection 
can be calculated, by using ( |47] >, as P X<» ( C H- Conse ‘ 
quently, the average probability of detection experienced on 
the fcth branch of the scenario tree in T stages is derived as 

(61) 

We aim at maximizing the average probability of detection 
calculated over all possible scenarios. Therefore, based on the 
scenario tree developed, we now represent our performance 
optimization in the following form 

JV„ 

ma ^r 5Z 7r(fe)P d|V> W (^ fc) ) (P2) 

s.t. 0 d ci k) d 1 Vfe, t\ 

! T • d fc) < (1 - vh Vfe, t; 

C ( 0 k) = C ( o j) Vk,r, 

C[ k:> = Ct ] whenever 
['0i = bPi 

where 7V st is the number of branches in the scenario tree. The 
last two constraints guarantee that the decision at each stage 
is non-anticipative , which means that, it does not depend on 
the observations made at subsequent stages. It is clear that, 
the decisions cannot depend on observations that are not yet 
available. 

As a significant advantage of © over ([Pi}, now we 
can decompose ( |P2| > into smaller (deterministic) optimization 
problems. Specifically, since the feth summand in the objective 
function of ( |P2| > only depends on C^ k \ we can convert ( |P2[ > 
into N st smaller optimizations, i.e., for k = 1,.... we have 

max7r (fc) P d | V , ( fc ) (V fc) ) CP3) 

s.t. 0 A d fe) =< 1 Vfe, t; 

! T • Ci fc) < (1 - rj)n Vfe, t; 

Co fc) =Co' ) Vfe, j; 

Ct^ = Ct ] whenever 

Note that ( |P2| ) is solved when all the N st optimizations in 
€1} are solved. Clearly, since iV sf grows exponentially with n, 
optimizing the system performance over the entire tree exhibits 
prohibitive complexity. However, decomposing the original 
optimization into a number of separate nonlinear programs 
enables us to propose a suboptimal approach for with 
much lower computational burden in terms of the number 
of nonlinear programs to be solved. Specifically, instead of 
considering all the N st optimizations, one may chose to solve 
only the problem which corresponds to the branch that has just 
been observed. In this manner, the decisions are made based 
on an average of the performance metric estimated by using 
the observations in T consecutive stages. 

Even dealing with one of the optimization problems in ( |P3| ) 
is computationally demanding in the sense that, its objective 


function is obtained by summing over nonconcave functions. 
We develop an affordable suboptimal solution by simplifying 
the optimization as a two-stage stochastic program, which 
will serve as the foundation for a very effective optimization 
approach developed in Section [V] Specifically, we truncate 
the scenario tree with T = 1 which leads to the following 
nonlinear optimization 

maxP*“L U{ k) ) (P4) 

£(fc) d|V>i V / 

s.t. 0 v Cl* 0 V 1 
1 T ■ c[ k) < (1 - I])n 


which is equivalent to 

Q-\a)^pl'E {HoM p L - 

^/pl'£ {Hl , hL} p L 

s.t. 0 d Pl d 1> l i ■ Pl < (1 - rj)n 


p L = argmin 

Pl 


(P5) 


Since the Q-function is monotonically decreasing, we have 
removed it form the objective function and turned the problem 
into a minimization. 

In the proposed cooperation scheme, for each different 
realization of the random interruptions we have a different set 
of cooperating nodes and consequently, we have a different 
detector. ( |P5| > represents a decision making process in which 
for each b/ the value of p/, is specified by applying the 
NPC on the detector corresponding to b /. Therefore, in case 
|P5} is solved for each b ^ generated by the random number 
generator, the Bernoulli distributions obtained are optimal for 
each realized cooperative detector in the NP sense. However, 
the decisions made by solving e are obtained based on 
only one set of observations. In Section [V] we discuss how 
to improve this decision making process by including more 
information about the effect of the random interruptions. 

( |P5| ) can be solved by applying the Karush-Kuhn-Tucker 
(KKT) conditions as follows. Note that the objective function 
in ( |P5| > is homogeneous of degree zero with respect to p [ J . We 
first relax the constraints and minimize the objective function 
to obtain a lower bound on the optimal value of ( |P5| ). Then, 
by using the homogeneity property of the objective function, 
we scale the derived point to find a feasible solution which 
attains this lower bound. More specifically, the lower bound is 
found by considering the Lagrange dual problem and applying 
the KKT conditions on the relaxed problem which yield [j 

C = S« 0 /2 [Q- 1 (a)I + KA]- 1 c (62) 

where Q is the global minimum of the relaxed problem, A = 
and c = E.j^ /2 q h ^, while n is the root of the 
polynomial equation [0 


[Q 1 (a)I + kA] X c 


= 1 


(63) 


and satisfies 

Q -1 (a)I + kA > 0 (64) 


2 The lower bound can also be obtained by the quadratic programming 
approach proposed in J7J. 
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Fig. 4. Scaling £ to a feasible point of Points #1 and #3 are mapped 
onto #2 and #4 respectively. The feasible point in each case is chosen to lie 
on the border of the feasibility region which is shown as a shaded area. 


Then, the optimal solution of (|P5|) is derived as 


interruptions. In consequence, by modifying ( |P5| ) we obtain the 
following optimization 

<3 _1 (“)\/ wrc fl|Wo w - a Tw 

P l = argmin- V , =- (P6) 

PL yw T c„| Wl w 

s.t. 0 A p L A 1, 1 T • PL < (1 - rj)n 

where 

a = m ~ t*u\H 0 = Z* Tp L ( E Kl^i] - E KM 

= Cf iU C“^P L (E[vl l \Hi] - E [u L \n 0 ]) 

» C^ U P| (E [ulI^x] - E [u L \Uo]) (68) 

and for h = 0,1, 


where 


Pl=AC 


A = 


r-C 




if A 1 

otherwise 


(65) 


( 66 ) 


Fig. HI demonstrates the scaling process indicated by ( |65j ) 
and d66f In this figure, optimal points #2 and #4 are derived 
by multiplying points #1 and #3 with their respective values 
of A. The feasible region of ( |P5| > is depicted as the shaded 
area. As seen in this figure, A in ( |66| maps £ to a point on 
the border of the feasibility region of ( |P5| >. This particular 
scaling is justified by recalling that, we solve as a fast 
suboptimal approach for optimization of a system whose best 
performance is established originally as Even though the 
optimal value of ( |P5| ) does not change by altering the norm of 
£, choosing a feasible point with the greatest norm leads to a 
better solution since the overall system performance depends 
on the norm of the decision vector p* L in general. 


V. Deflection Criterion 


Performance optimization based on deflection criterion 
means maximizing the DC which is defined as the variance- 
normalized distance between the centers of two conditional 
pdfs of the global test summary, i.e.. 


A KK»)lK]-iK(KK] 

v/Var {5(u)|7f 0 } 


(67) 


In this section, we show the link between the two-stage 
stochastic programming ( |P5| > and the DC. 

In the stochastic programming method developed in Section 
m the effect of the uncertainty is captured into the decision 
making process in the form of the statistics of the received 
reports conditioned on h L , i.e., bf} and C„| {Hfeibi} 

(for h = 0,1). We can improve this optimization approach by 
providing the decision maker with more information about the 
effect of 8 /. Specifically, we propose to make the decisions 
based on the statistics Pu\n h anc * C n |- Hfc , rather than their b/\- 
conditioned counterparts. Since the unconditioned statistics are 
the average of the conditional ones over Ol, they provide the 
decision maker with more information about the impact of the 


Hh ^yLu^yi ^yz. ^yi u 




= C K P L 


Pl 

(■-r.) 




94 


ulI Hh 


94 Ul + P L C Ul P h + cr?I 
hPrC 


-1 


U L | 'Hi 


P L+all 


94,, 


P lC Ul P l + I 


PlC 


pi f 

W L u r L u u L u 


L'^up.u 

(69) 


By using these approximations we see that the proposed 
optimization <[P6) is almost equivalent to 


s.t. 


P l= argmax = 

PL y / w T C fl |- Hl w 

0 A p L A 1, l 7 p L < (1 - r))n 


(P7) 


Since a 7 w = E [S(u)\Hi] — E [S(u)|H 0 ] an d = 

Var[S , (u)|'H/ l ], we clearly recognize the objective function in 
( |P7[ > as the DC. 

Hence, we have shown that modifying the decision making 
process in ( |P5| > to make the decisions based on overall (rather 
than conditional) statistics of the received interrupted reports 
is almost equivalent to optimizing the proposed system per¬ 
formance based on the deflection criterion. In addition, by the 
proposed modification the optimization variables derived are 
unpaired from 0 This means that unlike ( |P5| >, this problem 
needs to be solved only once. 

Now we proceed with developing a solution for m By 
using the approximations in ( |68j ) and ( |69| ). changing the 
variables, and with some further algebra we rewrite as 

77 ^ (P8) 


7r* = argmax ■ 


s.t. 


0 A 7r A D^ 2 !, 1 t D ct 1/,2 7t < (1 — rj)i 


where 7r = D<!/ 2 pL and D = 1 while and D CT are 

diagonal matrices whose elements are derived as 


= K 


(E [uM-EIvllIHo]) 1 


J 1,3 

(70) 

(71) 


(D a ) itj ± S itj (C UlU ww r C: [U ). . 

It is worth noting that, <(PHJ> is a convex-over-convex frac¬ 
tional program over a convex set. This type of optimization 
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is referred to as nonconcave fractional program and lacks 
a closed-form solution in general and is globally solved by 
using a special type of the BnB algorithm introduced in ]42| . 
Although the BnB is known to be very effective for small- and 
moderate-size problems, its high (worst-case) computational 
complexity, which grows exponentially with the number of 
optimization variables, might make it difficult to apply when 
a large-scale network is considered. 

In order to derive a computationally-affordable solution, we 
convert © into two nested optimization procedures. The 
inner procedure includes maximizing a QCQP and the outer 
procedure is realized as a one-dimensional optimization. Note 
that the feasible region of ( |P8| ) defines a polytope. The idea is 
to solve the problem on the intersection of the polytope with 
a ball centered at the origin, and to sweep the ball radius to 
scan the entire polytope searching for the global maximum. 

We formulate the QCQP by introducing a new constraint 
to make the feasible region of ( |P8| ) be on a ball of radius r, 
0 < r < y/Tr (Do-), and solve the following optimization for 
a constant r 

tp*(r) = -max7r T D7r (P9) 

r 7t 

s.t. 0 A 7r A Dy 2 l, l 2 1,/ “7 t < (1 — r])n, 

INI = r 

Then, the outer procedure, i.e., scanning the polytope, is 
realized as the following optimization 

r* = argmax ip* (r) (P10) 

r 

s.t. 0 < r < V /Tr (D CT ) 


Since ( |P10| ) is a one-dimensional problem, the computational 
complexity of the proposed nested loop is determined by §P9}. 
Solving © requires maximizing a convex quadratic objective 
function over a convex set, which is NP-hard in general. 
Fortunately, numerous relaxation techniques are available in 
the literature for solving nonconvex quadratic programs. A 
variety of these relaxations can be formulated as semidefinite 
programs (SDPs) and it is shown in [43] that the Shor relax¬ 
ation [44|j when strengthened with the so-called reformulation- 
linearization technique (RLT) |j45]| dominates other commonly- 
used semidefinite relaxations with optimality gap averaging 
3%. Accordingly, we first reformulate ( |P9| ) by introducing a 
new variable as V = 7T7r 2 which leads to 


maxTr(DV) (P9-a) 

V ,7T 

{ 0 A tt A Bl /2 ■ 1 

l r D c ( 1/2 7r < (1 - ri)n 
Tr(V) = r 2 
V = 7T7 r T 

The SDP relaxation is then realized by replacing the last 
constraint in ( |P9-a| > by V A 7T7 t t , while the RLT is con¬ 
structed based on using products of upper and lower bound 
constraints on the original variables to obtain valid linear 
inequality constraints on the new variable V. Consequently, 



Fig. 5. The proposed system performance in different SNR regimes and with 
different SNR dispersion levels among the sensing nodes. 


the resulting SDP is 

maxTr(DV) (P9-b) 

V, 7T 

l 2 1/,2 7r < (1 — rj)n 

V — 7rd T — d7r T + dd T > 0 
s.t. < V — 7rd r < 0 

V > 0 

V A 7T7T 2 

where d = D^ l. Note that the last constraint in this problem 
is equivalent to 


Therefore, ( |P9-b| ) is a convex optimization problem which can 
be solved, to any arbitrary accuracy, in a numerically-reliable 
and efficient fashion. 

This relaxation technique is known to provide very good 
approximate solutions for nonconvex quadratic programs in 
the form of ( |P9[ >. Nevertheless, ( |F9[ > can be solved for a global 
optimum in computational time complexity by a combinato¬ 
rial optimization algorithm developed in [46]. Therefore, the 
proposed approach in converting the fractional program ( |P8[ > 
into a set of QCQPs is a practically-appealing approach even 
if the global optimum is concerned. In the following section, 
the effectiveness of the proposed method in optimizing the 
proposed system performance is verified by simulation results. 

VI. Results 

In this section we visualize the performance of the proposed 
system by simulation results and make comparisons with the 
existing cooperative sensing schemes. In particular, we use the 
optimal linear combining scheme in as the benchmark. 

First, we consider three sensing nodes, i.e., I\ = 3, operat¬ 
ing at different SNR regimes and evaluate the detection and 
false alarm probabilities of the proposed system and compare 
them with the results predicted by our analysis. 

Fig. 0 depicts the detection probability under different 
average SNR levels and with different SNR dispersions among 
the sensing nodes. More specifically, the SNR levels consid¬ 
ered are SNRo — A, SNRo, SNRo + A where SNRo represents 
the average SNRs among the sensing nodes and A denotes 


7T 

V 


A 0 


(72) 
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Fig. 6. Effect of the CFAR threshold set by {46} on the false alarm rate of 
the proposed system. 


the SNR dispersion. The reporting channel noise variance is 
a 2 z = 10. The correlation coefficient p between the local 
sensing outcomes when Hi is true is 0.1 and the number of 
samples N used by each sensing node for energy detection 
is 20. Spatio-temporal correlations between the local sensing 
outcomes are realized by applying a moving average filter on 
the local sensing results which are realized under the AWGN 
with fixed channel gains corresponding to the different SNR 
levels. The Bernoulli distributions are considered as p/. = p () l 
where p (l indicates the average energy consumed for the local 
sensing and reporting purposes. For each set of theoretical 
and simulation results two values for SNR dispersion and two 
values for L are considered. Hence we have eight curves in 
Fig. 0 The results in this figure are obtained for significance 
level a = 0.01 and by 100,000 sample realizations for each 
average SNR level. 

It is clear from Fig. [5] that, regarding the probability of 
detection, the system performance predicted by the proposed 
analysis is in close agreement with the simulations. Moreover, 
we observe that the system performance is improved by taking 
into account the temporal correlations besides the spatial 
correlations. Specifically, the detection probability obtained in 
each case for L = 1 is higher than the performance obtained 
when L = 0. This is a reasonable observation as by using a 
more powerful compensation process, the degradations caused 
jointly by the interruption process and reporting channel 
contaminations are mitigated more effectively. 

Analytical and simulation results concerning the false alarm 
probability of the proposed cooperative sensing are depicted 
in Fig. [6] This figure depicts the effect of the threshold 
adaptation scheme in ( [46] ) on the false alarm probability (or 
more precisely, false alarm rate) observed. Specifically, in 
this figure, the false alarm rate corresponding to the adopted 
threshold level is compared against the false alarm rate that 
the system shows when using that threshold for detecting 
the PU signal. In this simulation, four cases are considered 
corresponding to po = 0.6, 0.8 and L = 0,1. We see in all 
these cases that, the proposed threshold adaptation works well, 
i.e., the observed false alarm rate of the detector is very close 
to the rate predicted by the analysis. 

In Fig. [7] the system performance under different levels 
of energy efficiency is compared with the performance of 


Fig. 7. CROC curves demonstrating the proposed system performance 
operating with different energy efficiency levels compared to the optimal linear 
fusion scheme. 

the optimal linear combining scheme developed in [6j—[8j. In 
this simulation, the complementary receiver operational char¬ 
acteristics (CROC) curves are plotted indicating the system 
performance in terms of the probability of false alarm and the 
probability of missed detection. The interrupted cooperation is 
simulated in two cases, i.e., for L = 0 and for L = 1, and in 
each case, four levels of energy efficiency are considered, i.e., 
efficiency levels of 30%, 50%, 70%, and 90% for three sensing 
nodes operating at SNR levels {12, 5,10} in dB. In this figure, 
we first observe the tradeoff between the energy efficiency 
and the detection quality. Specifically, the missed detection 
probability is increased when the system works with higher 
levels of energy efficiency. Moreover, for a given efficiency 
level, the proposed system exhibits a better performance for 
L = 1 compared to the case with L = 0. This better 
performance, which is also observed in Fig. [5] is justified by 
the fact that, a better compensation process leads to a higher 
detection probability. Moreover, the results in Fig. [7] indicate 
that, by using a more powerful compensation process the 
system performance becomes more robust to the degradations 
caused by the interruptions and reporting channel contamina¬ 
tions. For instance, we observe that the system performance 
degradation caused by increasing the energy efficiency level 
from 30% to 50% is smaller when a compensator with L = 1 
is used compared to the case in which a compensator with 
L = 0 is considered. This is the advantage of using a more 
powerful compensation process at the FC. 

Since in the optimal linear fusion scheme {6j-(8j only the 
spatial correlations are used at the FC, the CROC curves 
designated by L = 0 in Fig. [7] can be used to compare 
the proposed system performance with the performance of 
the optimal linear fusion. Specifically, it is clear in Fig. [7] 
that, the optimal linear combining performance is closely 
achieved by the proposed system. However, this performance 
is obtained while the energy consumed for cooperative sensing 
is significantly reduced by the proposed architecture. The 
energy efficiency level obtained is even more significant when 
temporal correlations are also taken into account as demon¬ 
strated by CROC curves corresponding to L = 1. 

In Fig. [8] the results of converting the fractional program 
m into a set of low-complexity SDPs is evaluated. We have 
considered three sets of nodes, operating at SNR levels as 
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Fig. 8. Performance of the proposed SDP for the interrupted cooperation 
with 30% energy efficiency compared with the performance of the optimal 
linear fusion scheme for K = 3, 9,15,30. 


Fig. 10. SNR loss vs. energy efficiency for Pfa = 0.1 obtained via the 
proposed cooperative sensing scheme with K = 3 sensing nodes operating at 
two different average SNRs, i.e., SNRo = 5,10 in dB, with dispersion levels 
as A = 0,1, 3 in dB. 



Fig. 9. Effect of imperfections in the available CSI on the probability of 
error in the proposed system and optimal linear fusion scheme. 


{12,5,7} in dB. The efficiency level considered is 30%. The 
CROC curves in this figure correspond to the cases where 
there are one, three, five, and ten nodes in each set of sensors. 
In this way, we investigate the system performance obtained 
by the proposed low-complexity suboptimal approach while 
increasing the number of sensing nodes. In addition, for each 
case the CROC curve corresponding to the optimal linear 
fusion is plotted for comparison. We can see in Fig. [8] that 
the results obtained by the proposed relaxation approach are 
very close to the overall detection performance obtained by 
the optimal linear fusion scheme. 

The effect of imperfect CSI on the overall detection per¬ 
formance is illustrated in Fig. [9] The proposed system perfor¬ 
mance is compared against the performance of the optimal 
linear combing scheme operating with the same erroneous 
CSI. The erroneous estimated first- and second-order channel 
statistics fi u and C u ^ h (for h = 0,1) are realized as 
normally-distributed random variables centered at their corre¬ 
sponding actual values while different estimation qualities are 
realized by changing the variances of these normal random 
variables. Three nodes operating at SNR levels {12,5,7} in 
dB are considered in this simulation. These nodes generate 
independent local sensing results as in 0, 0 - Assuming 
independent sensing outcomes leads to diagonal covariance 
matrices, i.e., less number of parameters to be estimated, 
and consequently, faster Monte Carlo simulations. Since im¬ 
perfections in the available CSI affects both the false alarm 


and detection probabilities, we compare the two cooperation 
methods by inspecting their probability of error defined as 
Pe = Pr{Ho}P/ + Pr{Pi}(l - P d ). In Fig. § P e is depicted 
vs. the normalized variance of estimation error, which is the 
variance of each parameter to be estimated devided by its 
squared value. The results show that as expected, performance 
of both detectors are degraded when the CSI error increases. 
However, the proposed interrupted cooperation outperforms 
the optimal linear fusion scheme when operating with im¬ 
perfect CSI and particularly when the error in CSI channel 
estimation is relatively high. 

The effect of SNR dispersion among the sensing nodes can 
be found in Fig. [TO] In this figure, the set of SNR levels experi¬ 
enced by the sensing nodes is {SNRo + A, SNRo, SNRo — A}, 
where SNRo denotes the average SNR of the sensing nodes 
and A represents the SNR dispersion. Two cases for the 
average SNR are considered, each with three values for the 
dispersion. That is, SNR 0 = 5,10 in dB and A = 0,1,3 in 
dB. It is clear that, a higher value for A represents a greater 
difference in the conditions in which the sensing nodes work. 
In order to compare the proposed system performance with 
the performance of the optimal linear combining, we consider 
the loss in SNR defined as 


SNR loss (dB) = SNR e (dB) - SNR W (dB) 


(73) 


where SNRe and SNR W denote the minimum SNR required 
at the SUs, incorporating the proposed method and optimal 
linear combining respectively, to meet a certain probability of 
false alarm, i.e., Pf = a, while having the same probability of 
missed detection. 

Fig. [TO] depicts the SNR loss vs. energy efficiency for 
a = 0.1. As shown in this figure, by increasing the efficiency 
level, the SNR loss is increased. This behavior illustrates 
the discussed tradeoff between the sensing quality and the 
energy consumption. In addition, we see that the SNR loss is 
decreased as the SNR dispersion increases. Therefore, a better 
performance is achieved when there is a greater difference 
in the conditions experienced by the SUs, i.e., when there 
is a higher need to discriminate between different sensing 


outcomes. In all cases depicted in Fig. 10 we observe that, 
significant energy efficiency can be achieved while the de- 
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tection performance is not heavily compromised. Moreover, 
the obtained energy efficiency is more significant when the 
dispersion is high. This is due to the fact that, when the SNR 
differences are high, the sensing outcomes from the nodes with 
low SNR do not considerably improve the overall detection, 
and therefore, the system suppresses their contribution by 
forcing them to go to the sleep mode more frequently as the 
efficiency level increases. Since the other nodes are relatively 
reliable, the overall detection performance is not affected 
significantly by this sleeping process. 

As another important observation, there are two breaking 
points on each curve in Fig. [TO] By increasing the efficiency 
level up to the first breaking point, the proposed system 
suppresses the adverse effect of the hidden node on the overall 
sensing process. Up to this point, the degradations caused 
by the interruptions and reporting channel are compensated 
by the linear estimator at the FC. The first breaking point 
is nearly where we have established an appropriate level of 
discrimination between the reliable and unreliable nodes and 
any further power reduction degrades the system performance 
since it affects the reliable nodes as well. Nevertheless, the 
compensator is still able to alleviate, to some extent, the 
degradations and consequently, the SNR loss does not grow 
as fast as it does when the second breaking point is passed. 
The second breaking point is nearly where the received 
reports quality is so low that the compensator is no longer 
able to reduce the performance degradation. Therefore, any 
further power reduction directly leads to significant loss of 
the detection quality. 


VII. Conclusion 

In this paper, a novel energy-efficient structure for spectrum 
sensing in CRNs has been proposed based on making random 
interruptions in the cooperation process among the CRs. The 
proposed system has been thoroughly modeled and analyzed, 
and an optimization problem has been developed in order 
to formulate a tradeoff taking into account the energy con¬ 
sumption at the local sensing and reporting processes jointly 
with the overall detection performance. Analytical solution of 
the optimization problem and the presented numerical results 
demonstrate that, significant levels of energy efficiency can be 
achieved by the proposed architecture. This energy efficiency 
is due to the fact that, unlike in existing cooperative sensing 
schemes, in the proposed design the discrimination between 
the reliable and unreliable nodes is obtained while no energy 
is wasted. Moreover, sensitivity of the overall detection to 
degradations in the local sensing and reporting processes is 
significantly reduced by the proposed architecture, leading to 
a more reliable cooperative spectrum sensing. 

Appendix I 

Second-Order Statistics of the Local Sensing 
Outcomes 

In the following, we derive C Ul and in terms of 

the listening-channel and PU-signal statistics. C lljU can be 
obtained similarly. We have 

(CuJij = ( R u L )„ - E [u k (m - l )] E [u n (m - r)] (74) 


where 


E [u k (m - 0] 
pv-i 

= E 


Li'=0 

N—l 


T | hks(N(m - l ) 

i'=0 

= yE\\h k s(N(m~l) 

o 

( E [g k ]\i 


i'= 0 

= N 


|s| 2 + tJ 2 


(75) 

i') + u k (N(m - l) - i') f 
i') + v k (N{m | 2 


and similarly, E[u n (m — r)\ = N yE[g n \ |s| 2 +cr 2 J, where 
g k — \hk\ 2 , gn — | h n | 2 are the channel power gains and 
|s| 2 = E[s 2 (m)\. Moreover, 


(R-uJij = E [u k (m - l)u n (m - r)] 

,2 


N -1 N-l 

= yyE \x k (N(m - l) - i')\ 2 I x n (N(m - r) - j' 

i'= 0 j'= 0 
N-l N-l 

= yyi E [9kg n ] R\s\* (N(l —r) + i'— /) 
i'—O j'= 0 

+ a 1 \ s \ 2 (E \g k \ + E [<?„]) + cr^ (l + S k=n 6 N (l-r),i'-j')} 


(76) 

2 


in which R\ s \ 2 (l — r) = E[\s(m — l)\ 2 \s(m — r)| 2 ]. We have 
assumed, without loss of generality that, the signals are wide- 
sense stationary (WSS). In addition, the noise variances are 
assumed to be the same for all sensors i.e., a Ui = cr v for all 
i. Note that 9f U| is a diagonal matrix whose main diagonal 
equals to the main diagonal of R Ul . 
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